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We develop the partitioning technique for quantum discrete systems. The graph consists of several 
subgraphs: a central graph and several branch graphs, with each branch graph being rooted by an 
individual node on the central one. We show that the effective Hamiltonian on the central graph 
can be constructed by adding additional potentials on the branch-root nodes, which generates the 
same result as does the the original Hamiltonian on the entire graph. Exactly solvable models are 
presented to demonstrate the main points of this paper. 
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I. INTRODUCTION 



The Schrodinger equation lies at the heart of quan- 
tum mechanics. Secular equation has analytic solutions 
only for a few very special cases. Approximation tech- 
niques and computational methods have been developed 
for treating such problem. Many of them are rooted in 
the partitioning technique [l^H which was introduced 
by Feshbach Q and Lowdin 4| independently. Discrete 
models, including quantum networks, have been a cor- 
nerstone of theoretical explorations due to their analyt- 
ical and numerical tractability [|[ , the availability of ex- 
act solutions, and the ability to capture counter-intuitive 
physical phenomena, such as non-spreading wavepacket 
Q and Bloch oscillation LZHp in linear chain. In recent 
years, optical lattice [Tol 111), photonic crystal [TH, [la |. 
etc. have increasingly permitted the experimental explo- 
ration of quantum discrete models. 

In this paper, we study the partitioning technique for 
quantum discrete systems. The concerned graph consists 
of several subgraphs: a central graph and several branch 
graphs, with each branch graph being rooted by an indi- 
vidual node on the central one. Applying the projection 
theory ji[ to such a graph, we show that the effective 
Hamiltonian on the central graph can be constructed by 
adding additional potentials on the branch-root nodes, 
which generates the same result as does the the original 
Hamiltonian on the entire graph. As the demonstration, 
we present two exactly solvable models, which correspond 
to real and imaginary potentials. 

This paper is organized as follows. Section [TTI shows a 
formalism for the partitioning technique in discrete quan- 
tum systems. Section Hill is the heart of this paper which 
presents a method to obtain the projection Hamiltonian. 
Section IIVI consists of two exactly solvable examples to 
illustrate our main idea. Section [V] is the summary and 
discussion. 
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FIG. 1. (Color online) (a) Schematic illustration of the graph 
consisted of two branch graphs a, b (cyan) and a center graph 
c (black). The dashed edges (red) represent the connections 
between them with A and B being the branch-root nodes, 
(b) The Lowdin's projection Hamiltonian for the center graph 
which is composed of the original Hamiltonian H c and addi- 
tional on-site potentials Va and Vb on the nodes A and B, 
respectively. 



II. PARTITIONING TECHNIQUE 

Lowdin has developed a partitioning technique in the 
algebra of matrices, with which various self-consistent 
field methods can be nicely formulated. In this proce- 
dure, the original Hamiltonian is simply transformed in 
a chosen discrete representation. The entire space is usu- 
ally divided into two subspaces, named a model space and 
an orthogonal space. The basic idea is to find an effec- 
tive Hamiltonian which acts only within the target model 
space but generates the same result as the original Hamil- 
tonian acting on the complete space [H,[H. The partition- 
ing technique enables us focus our interest on certain part 
of the system. In general, the effective Hamiltonian can- 
not be obtained explicitly, but provides a formalism to 
develop perturbation method. 

In the following we will show that, when the tech- 
nique is applied to a specific discrete system, the ef- 
fective Hamiltonian is of realistic significance. We con- 
sider a quantum graph, which is a collection of nodes 
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and edges. It is also equivalent to a single-particle tight- 
binding model. For simplicity, we partition the complete 
graph into three subgraphs, a central part c, two inde- 
pendent branches a and b. 

The Hamiltonian (or connectivity matrix) of such a 
graph has the form 



and more explicit form 



H 



H a H ac 
H ca H c H cb 
Hbc Ht, 



(1) 



where 



H b 

H c 



£ (4K) 6 01+h.c), 

£ (4K) c <j|+H.c.), 



(2) 
(3) 
(4) 



Here N, 



7— a,b,c 



denotes the dimension of the three sub- 



graphs. Kj~ a ' b ' c denotes the coupling between i and j 
of the graph 7, and reduces to the on-site potential for 
i = j. The connections between the subgraphs are 



N c N a 

i=l 3=1 
N c N b 

Hcb = Hi = -J2Y,ti\i)cb(j\- (6) 

i=l j=l 

where gjj~ a,h is the coupling strength between \j) and 
branch-root nodes \i) c . 

Our aim is the solution of the Schrodinger equation 



H\f k )=E k \f k ) 



(7) 



where 



1/*)= £ £/* 7 (oio. 

7— a,b,c /— 1 



(8) 



Then the Schrodinger equation can be written in the ma- 
trix form 



(9) 



H a 


H ac 




ft 


H ca 


H c 


H c b 







Hbc 


H b J 


u 




Halt H 


- H ac fk 


- E kfki 


(10) 






= E kfk> 


(11) 


Hbfi- 


I" H bc It 


= Ekfl 


(12) 



Under the condition of the existence of the inverse ma- 
trices (Ek — Hg)^ 1 and (E k — Hb)~, we have 



fk — (Ek - H a ) 1 H ac ft, 

fk = ( E k - H b ) 1 Hbcft, 



(13) 
(14) 



Then the Lowdin's projection Hamiltonian H c has the 
form 



where 



H c — H c + H a + Hb, 



H a — H ca {Ek — H a ) 1 H a 



H 



Hb — H c b (E k — Hb 



be ■ 



(15) 



(16) 
(17) 



Remarkably, the corresponding Schrodinger equation for 
the subgraph c (Eq. pTI) ) is reduced to 



Heft - E k ft, 



(18) 



i.e., formally H c can lead the same result as the original 
Hamiltonian acted with respect to the whole graph, then 
is referred as the effective Hamiltonian for central graph. 
Nevertheless, in general, one cannot treat Eq. (|18l) as 
usual since it is hard to obtain the explicit matrix form 
oiH r . 



III. EFFECTIVE HAMILTONIAN FOR 
CENTRAL GRAPH 

It can be seen from Eq. (fT5"|) that, H c is constructed 
based on the original subgraph H c . It indicates that the 
impact of two branch graphs can be projected on the tar- 
get graph as additional couplings or on-site potentials. In 
this paper, we investigate a graph with each independent 
branch graph connected to the central graph c via a single 
node on the central graph. This is crucial and our con- 
clusion is available for a graph with arbitrary branches. 
In the following we will show that H a and Hb have a 
concise form and clear physical meaning. 

The connections between the subgraphs are 

N a 

H ca =Hl = -J29*\A) ca (j\, (19) 
3 

N b 

Hob = Hl = -Y J 9 h j \B) cb {j\. (20) 
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Note that there is only one branch-root node for each 
branch, that is the unique restriction to the graph. 

We note from Eq. (fl9|) that the elements of H ca and 
H^ ac are all zeros except the row connecting to node A, 
i.e., 



H ca (to, n) = d mA 9ni H ac (m, n) = 6 nA (ff™)* ■ (21) 

Taking M a = (E k - H a )~ and assuming its existence 
for the considering eigenvalue Ek, we have 

H a (to, n) = E E Hca ( m >fi M ° UJ')Woc (j', n) 

3'=l 3=1 
N a N a 

j'=l 3=1 

N a 

= S mA S nA E 9]{9^)*M a (j, f) . (22) 

jd'=l 



Moreover, from Eqs. ([13]) and ([16)) we obtain 



Hcafk — H a fk 



(23) 



and its explicit form 



N a N a 

(J) = fk iA) E 9j(9f)*M a (j,f) ■ (24) 
i= i ij'=i 

Considering the non-trivial case (j4) ^ 0, the effective 
Hamiltonian H a can be expressed as 



H a (m, n) 



N a 



AOnA 



E^(j) 



(25) 



By a similar procedure we obtain expression for the ef- 
fective Hamiltonian Hf, 



H b (to, n) 



9-fk U) • 



(26) 



3 = 1 



Surprisingly, matrix H a (Hb) contains only one nonzero 
element H a (A, A) (Hi,(B, B)), which can be regarded as 
an effective on-site potential at the branch-root node A 
(B). Actually, this is caused by the unique restriction. 
Then the physics of the projection Hamiltonian is very 
clear: original target Hamiltonian with additional po- 
tentials at the joint sites. The effective potential is a 
weighted summation of the coupling strength {g]~°" b } 
and the corresponding amplitudes {fk(j)}- It would be 
noted that this conclusion can be generalized into graphs 
with more independent branches d, e, • • • . 



One can simply classify the branch graph as finite or 
infinite. For finite graph without flux, we have {g'j~ a ' b ' } 
and the corresponding {fk (j)} are all real, then the effec- 
tive on-site potentials are real. In contrary, for an infinite 
graph, when dealing with the scattering problem, the ef- 
fective on-site potentials could be complex. 



IV. ILLUSTRATIVE EXAMPLES 

In this section, two typical examples, which consist of 
finite and infinite branch graphs, are respectively inves- 
tigated to exemplify the formalism developed above. 

A. Finite chain 

We first take a finite chain N as an example, with the 
Hamiltonian in the form 



N—l 

ffchain = - J E d Z > (*+l|+H.C.). 

i=l 

It is well known that the eigenvalue Ek and the corre- 
sponding eigenvector fk are 



E k = -2Jcosfc, 
fk (j) = 
k = 



N + 
mr 

N+ 1 



1 



sin(fcj) . 



(27) 
(28) 



n e [1, N]. 



Now we divide the chain N as the central part N c and 
two branches N a , Nj, as mentioned above. The two 
branch-root nodes are located at the (N a + l)th and 
(N a + Ayth sites. From Eqs. (J23J) and ([25]). the pro- 
jection Hamiltonian can be obtained as 

JV„+JV C -1 

H c = -J E (NX* + l|+H.c.) (29) 

i=N a + l 

+V A \N a + 1) {N a + 1| + V B \N a + N c ) (N a + N c \ , 



where the on-site potentials are 



V A = 
V B = 



-J 



-J 



sin (kN a ) 
sin [k (N a + 1)] ' 
sin [k (N a + A^ e 
sin [k (N a + N c )} 



1)] 



(30) 
(31) 



In the Appendix IA 1\ it is shown that Ek is always the 
eigenvalue of H c and the corresponding eigenvector of T-L c 
accords with that of ifchain within the central chain c. It 
is noted that potential Va (Vb) does not exists in the case 
sin [k (N a + 1)] = (sin [ft (N a + N c )] = 0). Actually, the 
corresponding eigenfunction has vanishing amplitude at 
the node A (B), and Ek is also the eigenvalue of the 
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FIG. 2. (Color online) Schematic illustration of the concrete 
configuration for a scattering system. A ring as the scatter- 
ing center, connects to two semi-infinite chains L and R as 
waveguides with coupling —g. The wave function within the 
scattering center for a scattering state of the whole system is 
identical to an equal-energy eigen function of the projection 
Hamiltonian which is constructed by the center ring with ad- 
ditional on-site potentials Va and Vb added at the joint sites 
A and B. 



branch Hamiltonian H a (H ) simultaneously. From the 
viewpoint of the projection theory, the corresponding in- 
verse matrix (Ek — Ha)' 1 or (Ek — Hb)~ does not exist. 

Now we look at a concrete example in order to give 
a sense of the conclusion. Consider a 15-site chain 
with N a = 5, N c = 4, and N b = 6. Taking k = 
7r/4 as an example, the corresponding eigenvalue and 
eigenvector for the entire chain are E^/a = —*/2J, 
/tt/4 (J) = v2/4sin (jV/4), on the central chain c, 

= -(\/2, 1, 0, -l)/4. On the other hand, from Eqs. (|2?]>. 
([25]). ([3U]) and ([HI]), we have V A = - Jsin (5fc) / sin (6fc) 
= -V2/2J and V B = - Jsin (lOfc) / sin (9fc) = -y/2J. 
Then the corresponding effective Hamiltonian is 



B. Scattering problem 

In the above example, we can see that all the potentials 
are real. It was predicted that the infinite branches could 
induce the imaginary potentials. Here we are interested 
in the scattering solution of an infinite system. Quantum 
scattering and transport properties in quantum networks 
are important features in quantum information science 
d) EH EH • Now we consider an exactly solvable but non- 
trivial system to illustrate the main idea of this paper. 

The graph is constructed by a uniform ring system and 
two semi-infinite chains as the input and output leads, 
which is schemed in Fig. ([2]) . It is worthy to point that 
well-established Green function technique |15l - [l7| can be 
employed to obtain the reflection and transmission coeffi- 
cients for a given incoming plane wave. The correspond- 
ing wave function within the scattering center should be 
obtained via Bethe ansatz method. The Hamiltonian can 
be written as 

H s =H a +H b + H c (34) 
-V2J(|-l) QC (l| + |l) bc (AT + ll+H.c), 

with "Ha (Jib) represents a uniform input (output) waveg- 
uide as 



n a = ~Jj2 (l<-l).(i|+H.c), (35) 

i=-l 

+ 00 

^ = - J E(i j >^ + 1 i+ H - c -)> ( 36 ) 

i=l 

and the uniform ring as the scattering center is described 
as 



H f = - J H (I*) + + H - c - 16) (6I+V2 |9) (9|), 

(32) 

to solve H° s , we use the formula Eq. (|A6|) derived in the 
Appendix IA 11 It becomes 



sin (4k) (2 cos n —) = °> ( 33 ) 

which has the solutions E K = — 2 J cos k = V2J, 0, 
-V2J, and -3V^/2J. The corresponding eigenvec- 
tor for E K = —\[2J can be obtained as (/°/ 4 )^ oc 

(•\/2, 1, 0, — 1), which accords with wavef unction of whole 
system within the chain c, /°/ 4 . 



27V 

= -./^(|i) c (i + l| +H.c.) (37) 

i=l 

-V(\1) C (1\ + \N+1) C (N + 1\), 

where \2N + l) c = |l) c . 

There are on-site potentials V at the site |l) c and 
|JV + 1) , which are the two branch-root nodes, i.e., 
\A) c = |l) c and \B) c = |jV + l) c . The corresponding 
Lowdin's projection Hamiltonian depends on the energy 
Ek of the incident plane wave as well as the parameter 
V. To be concise, as an illustrative example, we would 
like to present the exactly solvable model, which is help- 
ful to demonstrate our main idea. Therefore, we will 
focus on the case: the incident wave has energy Ek = V 
£ (—2 J, 2 J). For such an incident plane wave, the scat- 
tering wave function can be obtained by the Bethe ansatz 
method. The wavefunction has the form, 
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/j?(0 = e tt < ,+1 >,Ze (-ao,-l], (38) 
fUl) = e m /V2, le [1,N + 1], (39) 
y*(I) = e «(i+^+i) j |g[l >00 ) j (40) 

where = /£ (2JV + 2 - Z). Then the effective Hamil- 
tonian "H a , can be obtained directly from Eqs. (|2"5)) 
and (|26p. which have the form 



Hamiltonian on the central graph can be constructed by 
adding additional potentials on the branch-root nodes, 
which can be expressed as a weighted summation of the 
corresponding wavefunction and generates the same re- 
sult as does the the original Hamiltonian on the entire 
graph. It indicates that the impact of the branch graph 
to the central one is local and takes the role of the on-site 
potential, A finite and an infinite exactly solvable models 
are presented to demonstrate our conclusion. 



H a {A, A) = - 
H b (B,B) 



V^Jfk ("1) 
V2Jf b k (1) _ 



= -2Je 
: -2Je lfe 



—ik 



(41) 
(42) 
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The projection Hamiltonian H c (H c = H c + 'Ha +Hb) is 



2N 



H c = -J t 

+2iJsink II) 



(43) 



[l\-2iJsmk\N + 1) C (N +1\ 



It is a VT symmetric non-Hermitian Hamiltonian. 
Since the seminal discovery by Bender fl8| . it is found 
that non-Hermitian Hamiltonian with simultaneous un- 
broken VT symmetry has an entirely real quantum me- 
chanical energy spectrum and has profound theoretical 
and methodological implications. In the Appendix IA 21 
it is shown the spectrum {s} of T-L c consists of a band 

Ej = -2Jcos(j7r/A), (44) 
( j G [1, N — 1] , 2-fold degeneracy) 

and two additional levels 



e± = ±V. 



(45) 



The eigenstates with eigenvalue ej can be decomposed 
into two kinds: symmetric and anti-symmetric with re- 
spect to the spatial reflection symmetry about the axis 
along the waveguides. For the scattering problem, only 
the symmetric states are involved. It shows that among 
the eigenvalues, the eigenvalue e + = V from the spec- 
trum {e} matches the energy E k (Ek = V) of the incident 
wave. Moreover, in the end of Appendix IA 21 it is shown 
that the corresponding eigenvetor for e + accords with /£. 
Thus it is in agreement with the conclusion of the parti- 
tioning technique that, there always exists a solution of 
the projection Hamiltonian to match the incident wave 
energy. 



V. SUMMARY 

In summary we apply the Lowdin's projection theory 
to the specified network, which consists of a central graph 
and several branch graphs. It is shown that the effective 



We acknowledge the support of the CNSF (Grant Nos. 
10874091 and 2006CB921205). 



Appendix A: Bethe ansatz solution 

In this Appendix, we will derive the central formula for 
studying the eigen problem of the projection Hamiltonian 
introduced in section HVl 



1. N-site uniform chain 

We consider a uniform chain with potentials at ends. 
The projection Hamiltonian is 

iv c -i 

H f = - J J2 + + H.c.)+Va |1) (1|+Vb |JV c ) (JV c 

(Al) 

where Va and Vb are defined in Eqs. ([50)) and (|3"T|) . The 
Bethe ansatz eigenvector has the form 



/,. ' • «,< " '-./ € [1,N C ] 



(A2) 



The Schrodinger equation H° s \f K ) = E K \f K ) can be 
written in the explicit form 



V A f K (1) - Jf K (2) = E K f K (1) 

-Jf* (J - 1) - Jf« (J + 1) = E K f K (j) 
3 e [2,N C -1] 
-Jf K (N c - 1) + V B f K (N c ) = E K f K {N c ) 



(A3) 



Substituting Eq. (|A2| into Eq. (fA3|) . we obtain 



J 2 sin [k (N c + 1)] + J {V a + V B ) sin {kN c ) (A4) 



-V A V B sm[ K (N e -1)] = 0, 

E k = -2Jcosk. (A5) 



Eq. (|A4[) determines the solution of k, while Eq. (|A5|) is 
the corresponding spectrum. Substituting Eqs. (|30[) and 
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(|3lp into Eq. (lA4j) . we have 



sin fc/ {sin [k (N a + 1)] sin [k (N a + N c )}} 
x {sin [k (N c - 1)] sin (nN c ) - sin (fcA c ) sin [« (A c - 1)]} 
+2 sin (kN c ) (cos k — cos fc) = 

(A6) 

which seems difficult to solve. However, it can be simply 
proved by straightforward algebra that, k = k is a so- 
lution for the equation. Accordingly, E K = — 2 J cos k = 
— 2 J cos A; is an eigenvalue of the effective Hamiltonian of 
Eq. (jAll) . Now we try to find the corresponding eigen- 
vector of E K . From Eq. (|A2[) . the first equation of Eq. 
[3| an d the expression of Va Eq. (f50|) , we obtain 



Bn 

A K 



-2ikN„ 



it indicates 



/ K (j)cxsin [k (N a+J )} 



which accords with the eigenfunction Eq. 
central chain N r . 



(A7) 



(A8) 



inside the 



f^U)={ ^ e i«(2iV+2-3) ^. _g e -t«(2iV+2-j) ' ( Al °) 

j G [N + 2, 2A] 

Substituting the above wave function into the following 
Schrodinger equation 



2* sinfc/ K (1) - f K (2) - f K (2 A) = EJ K (1) / J, 
-/« (j - 1) - /« (j + 1) = S re / K (j) / J, 

je [2,7V] U [TV + 2, 27V], (All) 
-/« (iV) - /« + 2) - 2t sinfe/ K (N + 1) 
= (N + 1) / J, 
after simplification, we obtain 



e *nN D _ e -iKN D 



A K e tK 
B K e~ iK 



= 0, (A12) 
£ k = -2Jcosk, (A13) 



2. Uniform ring as a scattering center 



where D± = sin fc ± sin k. 

The existence of the solution requires 



The projection Hamiltonian on a uniform ring is VT 
symmetric and can be expressed as 



2 A 



fff = -j£(|jMj + i|+H.c.) 

+2zJsinA:(|l) (1| - | A + 1) (A + 1|) , 



(A9) 



where \ j) — |2A + j). The parity operator V is given by 

7>|j) = ±|A + 2-j) 

and the time-reversal operator T obeys = —i. 

We note that the Hamiltonian also possesses the 
mirror symmetry with respect to the axis through the 1- 
th and (A + l)-th sites. This leads to the symmetric and 
antisymmetric solutions of the system. The symmetric 
Bethe ansatz eigenfunction f K has the form 



sin (kN) (sin 2 k - sin 2 k) = 0. (A14) 
The solution is 

K = nn/N,ne{l,N-l], (A15) 
K — k, 7r — k. 

the corresponding eigenvalue is Eqs. (|4*H 1431) . 

Obviously, E K = — 2Jcosk = — 2Jcosfc is an eigen- 
value of the effective Hamiltonian Eq. (IA9|) and the cor- 
responding eigenvecor is 



fn U) 



a ikj 



Jk{2N+2-j) 



je [l,A + i] 



, je [A + 2,2A] 



(A16) 



Therefore, the above eigenfunction f K accords with Eq. 
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